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ABSTRACT 


The  primary  purpose  of  this  report  is  to  present  analytical  tech¬ 
niques  for  the  prediction  of  instantaneous  parachute  temperatures  for 
any  prescribed  flight  condition  and  history  ot  operation.  Sample  calcu¬ 
lations  are  presented  which  illustrate  the  technique  by  comparing  measured 
parachute  temperature  histories  with  several  predictive  techniques.  A 
variety  of  exact,  transient  heat  conduction  solutions  are  included  (for 
a  cylindrical  geometry)  which  cover  a  range  of  trajectory  possibilities. 
This  was  done  because  the  transfer  of  heat  to  a  conventional  fabric,  or 
coated-fabric  parachute  is  governed  primarily  by  the  large  resistance 
to  heat  conduction  into  the  interior.  The  commonly  used  infinite  thermal 
conductivity  case  is  shown  r  >  produce  surface  temperatures  which  fall  con¬ 
siderably  below  the  surface  temperatures  predicted  by  more  exact  (complete) 
energy  balances.  For  the  extreme  applied  heat  flux  case  (atmospheric  en¬ 
try),  internal  temperatures  ar_  accurately  predicted  by  assuming  the  wall 
temperature  instantaneously  assumes  a  steady  state,  radiation  equilibrium 
value  (q^  *»  0)  and  internal  conduction  governs  the  process.  An  approxi¬ 
mate  method  is  presented  for  the  calculation  of  parachute  fabric  tempera¬ 
tures  when  both  the  internal  and  external  resistances  to  the  flow  of  heat 
are  of  importance.  The  results  are  compared  with  flight  test  data  on  a 
full-scale  parachute. 

The  distribution  of  this  Abstract  is  unlimited. 
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SECTION  J. 

it  is  well  known  that  by  «.Ite  aerodynamic  heating  effort  vehicles 
flying  at  supersonic  speeds  are  heated  to  such  an  extent  that  the  skin 
temperatures  aid  the  temperatures  inside  the  craft  require  the  special 
attention  of  the  designer.  In  addition,  a  growing  interest  in  aircraft 
vehicle  recovery  operations  has  directed  attention  to  a  variety  of  de¬ 
celeration  devices.  In  many  cases  a  sizeable  portion  of  the  initial  kinetic 
energy  of  the  primary  or  forsbody  is  converted  into  an  increase  in  in¬ 
ternal  energy  of  the  retardation  device — and  the  aerodynamic  heat  trans¬ 
fer  of  the  deceleration  device  becomes  a  design  factor. 

Aerodynamic  decelerates  must  possess  a  large  drag-to-weight  ratio 
and  normally  must  be  capable  of  being  stored  in  a  small  volume.  There¬ 
fore,  most  drag-producing  devices  are  relatively  thin  and  flexible.  Towed 
decelerators  such  as  ribbon  or  hyperflo  parachutes,  balloons,  halluces, 
etc.,  have  flexible  caonpies.  Other  proposals  incorporate  structural 
members  (cones  or  paragliders)  or  Hybrid  systems  such  as  the  AVEO  drag 
brake  (1)*  and  the  NASA  rotornel  (2). 

Impermeable  retardation  devices  are  aerodynamical ly  unstable  when 
utilized  in  a  supersonic  flow.  As  a  result,  nearly  all  the  retardation 
devices  discussed  above  employ  extended  porous  surfaces  n  conjunction 
with  other  venting  techniques  (see  Figure  1).  Most  conventional  textile 
libers  either  melt  or  are  serious  ly  degraded  by  exposure  to  temperatures  ab  .»ve 
It  is  desirable  to  rep  la  e  these  libers  with  materials  of  higher 

♦Numbers  enclosed  by  parentheses  refer  to  references. 


thermal  durability  such  as  fine  metallic  wires  in  which  operational 
temperatures  are  close  to  2000°F. 

It  i.o  felt,  that  for  the  purpose  of  analysis  the  flow  field  associated 
with  a  porous  retardation  device  Includes  large-scale  phenomena  enveloping 
the  entire  device  and  small-scale  processes  associated  with  the  local  flow 
about  a  single  element  in  the  mesh  material  itself.  The  large-scale 
phenomena,  such  as  parachute  shape  or  size  and  the  velocity  and  altitude, 
involve  the  external  aerodynamics  of  the  problem  and  contribute  the  end 
or  boundary  conditions  (such  as  pressure  ratio  or  Reynolds  number)  to  the 
small-scale  phenomena. 

The  mesh  element  geometry  under  consideration  is  a  porous  weave 
consisting  of  individual  round  wire  filaments.  Wind  tunnel  (3)  and  free- 
flight  tests  (4)  of  ribbon  parachutes  have  shown  considerable  degradation 
of  performance  due  to  the  high  heat  tran  ler  coefficients  near  the  edges. 
Lately  wind  tunnel  tests  or  small  models  have  indicate!  that  reasonable 
inflation  and  attitude  stability  may  be  attained  with  a  fir.e-meeh 
roof  material  (3).  The  flow  through  a  single  pore  is  taken  to  be 
analogous  to  the  flow  through  a  nozzle.  At  supersonic  flight  speeds 
and  at  very  high  porosities*  discrete  shock  waves  form  ahead  of  each 
cylinder.  If  the  porosity  is  lower,  i.e.,  insufficient  to  permit  super¬ 
sonic  flow  to  be  "started"  in  each  nozzle  In  a  manner  similar  to  the 

*Poroslty  is  defined  geometrically  as  the  ratio  of  open  area  to 
total  area  of  the  forward-facing  surface  of  the  retardation  device. 
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starting  process  in  a  closed-channel  wind  tunnel,  a  continuous  shock 
appears  ahead  of  the  entire  mesh  surface.  At  such  intermediate  poros¬ 
ities  and  sufficiently  large  Reynolds  numbers,  a  discrete  boundary  layer 
will  form  on  each  cylinder.  Only  at.  very  low  porosities  and  Reynolds 
numbers  will  the  individual  boundary  layers  merge  to  completely  fill  the 
pores  and  establish  a  continuous  boundary  layer  over  the  entire  mesh  sur¬ 
face.  Only  the  intermediate  porosity-discrete  boundary  layer  case  will  be 
studied  in  this  report. 

Large-Scale  Phenomena 

The  range  of  interest  in  che  present  study  consists  of  altitudes  from 
sea  level  up  to  200,000  feet  and  flight  Mach  numbers  ranging  from  3  to  6. 

A  typical  "problem"  might  be  Mach  5  operation  above  70,000  feet  with  a 
nominal  parachute  diameter  of  2  to  10  feet  and  geometric  porosities  of  25%. 

A  concave  Hemisphere  canopy  is  a  typical  shape  considered  for  supersonic 
operations  although  specification  of  the  exact  geometry  is  not  necessary. 

In  previous  work  \6) ,  it  was  felt  that  a  good  starting  point  for  study  of 
the  large-scale  phenomena  was  the  simple  concave  hemisphere.  The  experi¬ 
mental  results  of  Reference  6  indicate  that  at  angular  locations  up 
to  75  degrees  from  the  stagnation  point,  the  static  pressure  on  the  upstream 
surface  of  the  concave  hemisphere  is  equal  to  the  total  pressure.  This 


means  that  there  Is  practically  no  flow  near  and  parallel  to  the  surface 
of  the  hemisphere  e  cept  near  the  t  dge.  Of  course,  in  supersonic  flow,  a 


nearly  normal  shock  viould  stand  ahead  of  the  hemisphere  and  this  would 
produce  a  completely  subsonic  flow  regime  inside  of  the  concave  hemisphere. 
In  the  absence  of  any  parachute  porosity,  the  pressure  on  the  upstream 
side  of  the  chute  is  closely  approximated  by  the  total  pressure  behind  the 
normal  shock.  In  determining  the  pressure  on  the  rearward  (downstream)  face 
of  the  parachute,  we  are  led  to  a  base  pressure  phenomenon  described  by 
H.  H.  Korst  (7). 

The  discussion  in  the  previous  paragraph  illustrates  a  technique  for 
computing  the  pressures  and  hence  the  pressure  ratio  across  a  chute  with 
no  porosity.  The  analysis  for  the  base  pressure  problem  has  been  extended 
to  the  case  of  flow  through  the  chute  with  small  momentum  only.  Since 
this  last  condition  is  rarely  met  in  »_tie  case  of  practical  parachute 
configurations ,  che  pressure  ratio  across  the  chute  openings  must  be 
found  from  experiment. 


Small-Scale  Phenomena 

A  high- temperature,  high-density  flow  is  produced  inside  the  canopy. 
Tliis  flov»  passes  through  the  individual  openings  of  the  mesh.  Heat  is 
transferred  from  the  hct  gases  to  the  mesh  elements.  Therefore,  it  is 


•f- 

Schlieren  photographs  of  supersonic  wind  tunnel  tests  Indicate 
that  trie  normal  shock  one  anticipates  ahead  of  a  low-porosity  blunt  cano¬ 
py  Is  altered  into  a  conical  shape  bj  the  viscous  forebody  wake.  Mainly 
the  low  energy  portions  pass  through  the  canopy. 


the  distribution  of  heat  flux  to  the  surfaces  of  a  crossed  ii.esh  which 
is  of  primary  Importance.  The  typical  element  may  be  considered  as  a 
cylinder  aligned  normal  to  the  flow.  The  element  is  bounded  on  each  side 
by  an  opening  (slot)  through  which  the  oncoming  flow  passes.  The  ap¬ 
proaching  flow  may  or  may  not  separate  upstream  of  the  cylinder.  This 
process  has  not  been  clearly  defined. 

The  flow  passing  through  the  slots  exhibits  properties  which  depend 
upon  the  overall  applied  pressure  ratio.  In  most  flight  applications, 
the  applied  pressure  ratio  will  be  such  that  the  flow  in  the  slot  may 
be  sonic.  Additional  energy  is  available  for  further  expansion  of  the 
stream  as  it  emerges  from  the  slot  such  that  localized  regions  of  super¬ 
sonic  flow  will  occur  downstream  of  the  slot.  Between  the  supersonic 
jets  emerging  from  the  slot  and  the  rearward-facing  surface  of  the  cylinder, 
two  regions  are  found.  The  first  is  a  conventional  free-shear  layer  in 
which  the  streamwise  component  of  velocity  diminishes  from  the  value  found 
in  the  jet  to  a  near-zero  condition.  Adjacent  to  the  shear  layer  there 
exists  a  small  region  of  reverse  flow  with  a  circulating  vortex.  The 
local  pressure  In  the  separation  bubble  is  determined  by  local  base  pres¬ 
sure  phenomena  (the  flow  issuing  from  the  slot) ,  The  mean  near-wake  pres¬ 
sure  which  exists  further  donwstream  is  determined  by  the  gross  geometry 
and  the  freestream  Mach  number.  These  two  pressures  are  generally  not 
equal.  Therefore,  further  recompression  exists  downstream  of  the  jets  is¬ 
suing  from  the  Individual  slots.  The  local  base  pressure  behind  an  indi¬ 
vidual  ribbon  is  therefore  determined  by  the  slot  pressure  ratio  and  the 
recompression  process  in  the  jets. 
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Flow  Field  Analysis 

On  the  basis  of  the  preceding  discussion,  the  following  assumptions 
v/ere  made  regarding  the  flow  field  about  the  parachute: 

1)  No  consideration  is  given  to  the  presence  of  the  forebody  wake  and 

to  corresponding  t.  edifications  of  the  canopy  shock  wave.  The  local  flow 
parameters  (M^,  p^,  ,  p^)  ahead  of  the  body  (station  1)  are  determined 

solely  b2r  the  trajectory. 

2)  A  detached  shock  stands  in  front  of  the  inlet  area  (see  Figure  2). 
The  shock  is  assumed  to  be  a  normal  shock.  The  shock  is  not  influenced 
by  the  riser  lines.  The  fluid  properties  downstream  of  the  shock  (sta¬ 
tion  2)  are  determined  by  the  normal  shock  relations. 

3)  The  flow  through  the  skirt  is  one-dimensional.  The  local  flow  param¬ 
eters  are  determined  using  a  steady,  one-dii  -i  onal,  isentropic,  com¬ 
pressible  flow. 

4)  The  local  flow  through  the  roof  is  analogous  to  the  flow  in  the  sonic 
throats  of  a  multinozzle  grid.  The  sonic  pressure  in  the  mesh  minimum 
area  is  found  from  the  isentropic  expansion  from  the  fictitious  reservoir 
(Figure  3)  with  conditions  (Pq>  Tq)  that  exist  behind  the  shock.  Con¬ 
tinuity  of  mass  flow  through  the  inlet  and  throat  can  be  used  to  deter¬ 
mine  .e  approach  flow  velocity  inside  the  roof  at  station  3. 

5)  This  approach  to  calculating  the  heat  flux  rates  to  porous  retar¬ 
dation  devices  is  also  subject  to  the  following  limitations: 

a)  porosities  such  that  individual  boundary  layers  exist  on  each 
filament  and 

b)  continuum  boundary  layer  concepts  on  each  filament. 
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Figure  2.  Parachute  Flow  Field 
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Figure  Trajectory  Nomenclature 


Heat  Transfer  Phenomena 


Most  parachutes  are  constructed  of  fabrics  which  are  very  poor 
thermal  conductors.  Examination  of  parachute  canopies  after  experimental 
tests  has  revealed  the  Importance  of  parachute  material  mass  and  its 
distr'* huMor> .  Are-;  remote  overlapped  seams  have  obtained  higher 

temperatures  than  in  the  higher  mass  seam  areas.  The  thermal  heat  sir.'/ 
of  a  giver,  component  is  significant  in  determining  the  final  temperature 
attained  by  that  component. 

When  the  decelerator  is  being  heated,  the  maximum  temperature  the 
surface  can  attain  is  the  equilibrium  temperature — i.e.,  the  temperature 
resulting  from  a  balance  between  convective  heating  and  radiation  cooling. 
Therefore  the  maximum  temperature  which  may  develop  is  strongly  dependent 
on  the  film  heat  transfer  coefficient  ai.d  the  emittance  of  the  material. 

The  radiation  equilibrium  temperature  closely  approaches  the  adiabatic 

wall  temperature  at  low  Mach  numbers  where  the  total  temperature  is  snail. 

4 

The  departure  is  significant  above  Mach  number  2.0  due  to  the  T  variation. 

For  the  '‘actual"  or  "instantaneous"  temperature, the  local  heat  transfer 
coefficient,  the  surface  emittance,  and  the  thermal  diffusivity  of  the 
decelerator  material  must  be  known  for  each  point  of  the  trajectory.  Fre¬ 
quently  time-temperature  histories  are  computed  using  the  assumption  of 
infinite  object  thermal  conductivity  coupled  to  the  neat  sink  ability  of 
the  parachute.  Calculated  maximum  temperatures  involve  a  balance  between 
convective  heat  input,  and  radiative  he  a'  loss  and  do  not  include  heat 
storage  terms.  The  zeto  convective  heat  flux  wall  temperature — referred 
to  as  t lie  recovery  or  adiabatic  wall  temperature,  1 


results  in  wall 


temperatures  which  are  close  to  the  freestream  total  temperatures  since 

the  expression  for  T  is 

aw 


T 

aw 


CO 


1  +  r  (- 


1  +  (“ 


and  the  values  of  r  fall  in  the  range 


0.85 


(Pr> 


1/2 


r,  <  r  s  r 
lam  turb 


(Pr) 


1/3 


0.90 


(1) 


(2) 


Some  authors  have  recommended  porous  drag  devices  because  of  the  ad¬ 
vantages  offered  from  the  points  of  view  of  flexibility,  packaging,  manu- 

x 

* 

facture,  and  aerodynamic  stability.  There  is  a  serious  disadvantage  of 
increased  heat  transfer  from  the  surrounding  air  to  the  porous  retarda¬ 
tion  device.  When  the  retardation  device  is  made  up  of  a  porous  weave 
consisting  of  round  filaments,  and  the  poroc'ty  of  the  weave  is  such  that 
a  sing.e  normal  shock  forms  ahead  of  the  body,  the  total  heat  transfer  to 
such  a  woven  cloth  is  equal  to  the  sum  of  the  heat  transfers  to  the  indi¬ 
vidual  filaments.  It  is  a  well-established  fact  that  the  stagnation  point 
heating  rates  vary  as  the  inverse  square  root  of  the  element  radius.  There¬ 
fore  the  ratio  of  the  heating  rates  of  porous  materials  to  solid  materials 
is  equal  to  the  square  root  of  the  nose  radius  to  filament  radius.  For 
example,  a  weave  consisting  of  0.010-inch-diameter  filaments-50  to  the 
inch  (25  , crcent  geometric  porosity),  receives  50  times  as  much  heat  pur 
unit  area  as  a  solid  material  of  1  foot  radius  under  the  same  approacn 
flow  conditions.  The  area  available  for  heat  transfer  is  also  increased. 
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The  equilibrium  temperature  is  affected  to  a  lesser  amount  because 
woven  materials  have  higher  emissivities  than  solid  surfaces  of  the  same 
parent  material — a  geometrical  result  due  to  the  great  number  of  cavities. 
These  cavities  tend  to  alter  the  reflection  characteristics  and  therfore 
to  increase  the  absorptivity  above  that  of  a  solid  surface.  At  equilib 
rium,  the  emissivity  and  absorptivity  are  equal  and  incrc  ?sed  emissivity 
results. 

Mass  Transfer  Cooling  Effects 

At  the  present  time  useful  metallic  and  cevaai c  fibers  are  still  in 
the  developmental  stage  and  parachute  materials  can  only  be  used  below 
their  softening  temperature  (for  nylon  350°F)-  This  limitation  suggests  the 
use  of  non-permanent  coatings  that  sublime  or  ablate.  In  consideration 
of  high  -  temperature  materials  for  parachute  application,  coated  f  abrics 
or  coated  fibers  are  useful.  The  short  aerodynamic  heating  periods 
associated  with  first-stage  deceierators  permit  the  use  of  ablative-type 
insulating  coatings  of  conventional  internal  fibers  (nylon  and  Nomex) .  For 
orbital  re-entry  trajectories  coated  metal  fiber  fabrics  would  be 
effective.  It  has  been  demonstrated  (8)  that  subliming  coatings  do 
substantially  increase  the  operating  lifetime  of  the  base  material.  The 
increased  thermal  protection  results  primarily  from  the  ac.ded  mass  (heat, 
sink  effect).  The  extreme  mechanical  environment  (parachute  opening 
shock,  flutter,  etc.)  Is  extremely  detrimental  to  the  bond  between  the 
coating  and  parent  material  such  that  frequently  the  coating  material  is 


removed  hv  mechanical  rather  than  thermal  causes. 


Trajectory  He at  Transfer  Considerations 

The  greatest  source  of  uncertainty  in  engineering  calculations  of 
the  energy  exchange  process  to  the  mesh  elements  lieu  in  the  determina¬ 
tion  of  the  convective  heat  flux  as  characterized  by  a  heat  Lransfer 
parameter  Nu  which  is  a  function  of  a  number  of  parameters.  This  is 
expressed  by  the  following  dimensionless  relationship 

Mu  =  f (Re,  Mj,  Kn,  T^/'Tq,  geometry).  (3) 

The  Reynolds  number  Re  describing  the  flow  ccauitions  around  the 

fibers  will,  in  this  report,  be  based  on  cylinder  diameter  D  and  on  the 

average  velocity  and  properties  in  the  flow  through  the  slots  between 

the  mesh  elements  where,  for  supercritical  pressure  ratio,  sonic  velo- 

o 

city  exists.  This  Reynolds  number  is  denoted  as  Re*  =  r ; —  .  In 

u* 

Figure  4,  Ke*/D  is  plotted  versus  descent  Mach  number  for  altitudes 
ranging  from  sea  level  to  200,000  feet.  This  i  Reynolds  number  is  pre 
seated  since  it  does  not  require  th"'  specification  of  a  body  dimension 
which  varies  with  a  specific  parachute.  The  large  values  of  Re*/D 
which  ex'  low  altitudes  produce  high  heat  flux  rates.  The  short- 

dashed  curve  is  a  curve  of  constant  dynamic  pressure  (connected  with 
the  descent  velocity!  of  1  lb.  per  square  foot — below’  which  reliable 
sell-inflation  oi  flexible  drag  devices  cannot  be  expected.  The  curve 
labeled  "stagnation  temperature  =  1000°F"  illustrates  tiie  heat  transfer 
problem  as  determined  by  Mach  number  effects.  The  operating  regime  of 
conventional  parachutes  lies  within  the  envelope  of  the  two  preceding 
curves . 
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eynolds  dumber  Versus  Flight  ?'ach  dumber  for 
Various  Flight  Altitudes 


The  descent  Mach  number  has  only  an  Indirect  influence  upon  the 
heat  transfer  to  the  mesh  elements  insofar  as  it  determines  the  pressure 
and  teoiper'‘ture  both  inside  the  canopy  and  behxnd  it.  For  the  flow  and 
heat  transfer  around  the  ribbons,  the  pressure  rat-o  p^/p^  *8  a  moce 
suitable  parameter.  Its  use  has  the  further  advantage  that  it  can  be 
readily  measured  in  a  wind  tunnel  experiment  on  a  parachute  model. 

The  pressure  ratio  curve  presented  in  Figure  4  is  the  ratio  of  the  total 
pressure  behind  a  normal  shock  to  the  freestream  static  pressure  and 
represents  the  minimum  pressure  ratio  which  should  be  studied  at  any 
simulated  flight  Mach  number. 

The  Frandtl  number  Pr,  whose  value  is  approximately  0.7  under  normal 
atmospheric  conditions,  deviates  up  to  +  10  percent  from  this  value  up 
to  temperatures  of  approximately  6,000°R.  At  altitudes  up  to  400,000  feet 
the  stagnation  temperatures  for  Ma  £  10  will  be  below  6,000°R.  The  de¬ 
viation  of  the  Prandtl  number  from  the  value  0.7  will  be  neglected. 

Ine  Knudsan  number  Kn  wili  enter  as  a  parameter  only  at  very  low 

densities.  It  can  be  expressed  as  a  relation  between  Mach  number  and 

M1 

Reynolds  number,  for  instance  by  — —  .  In  this  expression  the  mean  free 

/Re 

path  of  the  molecules  is  related  to  the  boundary  ayer  thickness.  The 

M1 

determination  cf  the  critical  Value  of  — ™  depend-  upon  the  shape  of 

/Re 

the  body  which  affects  -he  boundary  layer  g/o'th.  However,  for  an  ap¬ 
proximate  evaluation  it  seems  justified  to  introduce  M*  and  Re*  as  charac¬ 


teristic  vaxues  in 


the  Knudsen  number.  Hermann  (9) 


give  as  critical  value 


10  .  in  the  present  case,  M*  =  1.  Thus, 


;ont Lnuura  flow  can  be  expected  to  exist  for  descent  conditions  at  which 


Re*  >  1000.  The  conditions  for  which  this  is  the  case  are  evident  from 


Figure  4.  The  results  presented  in  thj s  rep ort  are  valid  for  the  continuum 
fJow  regime.  Summarising  the  foregoing  discussion  on  the  Influence  of 
the  various  parameters  upon  the  heat  transfer  on  a  parachute  fiber,  we 
can  expecc  a  relation 


Nu  =  f (Re*, 


(4) 


This  relation,  then,  has  to  be  determined  by  experiments. 

As  mentioned  above,  it  is  convenient  to  defir, e  a  Reyno]  -"s  number 
based  on  the  velocity  and  flow  properties  in  the  sonic  'rifice  and  the 
slot  width 


Re*, 


p*a* 


(5) 


In  the  flight  case  the  static  temperature  T,  and  pressure  p,  are  deter- 

-X. 

mined  by  the  flight  altitude.  A  total  p"essure  pQ  and  temperature  TQ 

are  determined  once  the  flight  Maca  numoer  As  prescribed.  Since  the  bow 

shock  is  assumed  to  be  normal,  the  total  pressure  behind  the  normal 

shock  p*  (inside  the  canopy)  is  determined  from  pn  ai  1  M. .  The  total 

density  behind  the  bow  shock,  is  computed  from  p f . ,  T  and  the  proper 

u  U 

equation  of  st3te.  Finally,  the  sonic  density  and  speed  are  determined 
from  an  isentropic  expansion  from  p'q»  to  P*  a!K*  T*.  The  viscosity 

u*  is  a  known  function  of  tne  temperature  T*. 

Once  a  method  for  calculating  t h heat  transfer  coefficient  is 
selected,  an  appropriate  heat  balance  solution  must  be  applied.  Either 
or  two  approximate  proceuures  car.  be  used,  as  follows: 


1 .  One-dimensional  transient  t  glut  ion : 

Consider  the  fabric  wall  as  a  cylinder  exposed  to  aerodynamic 
heating  on  the  surface  A  heat  balance  then  can  be  written 
for  the  conditions  at  :he  outer  surface.  The  condition  for 
heat  input  at  the  outer  surface  by  convection  alone  from  a 


fluid  at  an  adiabatic  wall  temperature  (T  )  is 


h(T)  [T  -  T(r,t)] 

si W 


(6) 


This  equation  implies  that  the  element  is  subjected  to  a  uni¬ 
form  thermal  environment  around  its  periphery.  Thus,  one 
can  approximate  the  heating  of  the  roof  panel  element  by  a 
cylinder  subject  to  a  constant  temperature  environment  at  the 
instant  of  time  under  consideration.  As  a  first  approxima¬ 
tion,  the  radiation  away  from  the  outer  surface  is  neglected 
but  can  be  included  if  necessary  as  surface  temperatures  reach 
values  where  they  become  an  important  factor.  The  rate  of 
heat  conduction  into  the  cylinder  then  can  be  calculated  in 
conjunction  with  one  of  the  suggested  methods  for  calculating 
the  heat  transfer  coefficient,  by  iteration  over  a  finite  time 
interval.  This  equation  is  solved  over  the  critical  decelera¬ 
tion  time  with  aid  of  a  computer,  or  it  is  estimated  by  using 
thermal-response  charts  such  as  those  available  in  Section  IV, 
Case  C.  This  transient  type  of  heat  balance  solution  is  use¬ 
ful  particularly  in  situations  where  the  deceleration  takes 
place  over  a  short  time  interval. 
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2. 


Heat  balance  solution: 


Consider  the  fabric  wall  as  a  cylinder  exposed  to  aerodynamic 
heating  or.  its  external  surface  which  behaves  like  a  heat  sxnk 
with  practically  no  temperature  gradient  through  its  thickness. 
Then,  the  following  heat  balance  can  be  written: 


hS, 


(T  -  T  )  -  oS,eT 
aw  w  3  w 


p  V  c 
m  m  p 


m 


dT 

dx 


(7) 


This  heat  balance,  developed  in  Section  I,  is  written  at  the 
outer  surface  on  the  assumption  that  there  is  no  internal  ra¬ 
diation-energy  interchange.  This  equation,  combined  with  the 
method  selected  for  calculating  the  heat  transfer  coefficient, 
can  then  be  solved  by  iteration,  by  assuming  various  wall 
temperatures  for  the  instant  of  time  under  consideration  until 
a  balance  is  achieved  (i.e.,  a  quasi-steady-state  condition 
is  reached  along  the  trajectory  path  for  a  finite  time  interval). 
This  heat  balance  solution  is  useful  particularly  in  situations 
where  the  material  has  practically  no  temperature  gradient  with¬ 
in  it  (Newtonian  heating).  Thus,  a  simple  estimate  of  the  sur¬ 
face-temperature  rise  can  be  made  if  one  assumes  that  such  a 
temperature  gradient  does  not  exist. 

Either  Equation  (6)  or  (7)  may  be  used  to  estimate  the  maximum 
temperature  the  parachute  may  attain  during  the  heating  process.  This 
equilibrium  temperature  will  normally  exceed  the  actual  surface  tempera¬ 
ture  during  the  parachute  heating  cycle.  There  are  cases,  such  as  atmos- 


pheric  entry,  in  which  the  instantaneous  equilibrium  temperature  looks 
like  one-half  a  cycle  of  a  sine  wave.  The  complete  solution  of  Equation 
(7)  which  incorporates  the  internal  heat  sink  effect,  leads  to  actual 
wall  temperatures  which,  at  times,  exceed  the  equilibrium  temperature  as 
sketched  below: 


TIME  r 


The  methods  described  above  must  be  modified  if  the  decelerator 
reaches  its  charring  temperature  or  melting  temperature  before  attaining 
the  equilibrium  temperature.  Certainly  computed  temperatures  that  are  in 
excess  of  the  melting  temperatures  have  no  physical  significance.  When 
the  surface  temperature  reaches  the  melting  temperature  the  problem 
shifts  to  a  transient  heat  conduction  problem  with  a  constant  surface 
tempc i  iture. 
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SECTION  II.  EXTERNAL  ENERGY  TRANSFER 


Energy  Balance 

The  time- temperature  history  of  a  deceleration  device  is  governed 
by  the  principle  of  energy  conservation  which  is  expressed  by  the  relation 
between  energy  flow  rates  Q  as 

^entering  ^leaving  ^stored 

In  the  case  of  a  deceleration  device  moving  through  the  atmosphere  the 
foil  "wing  energy  flows  appear: 

1)  Energy  flow  into  the  parachute  by  convection, 

2)  Energy  gain  due  to  incident  solar  radiation,  Q 

so 

3)  Energy  loss  from  the  parachute  by  radiation  to  space,  Q 

rsp 

4)  Energy  loss  due  to  radiation  to  the  earth,  Q 

.  e 

5)  Energy  stored  by  the  parachute,  Q 

These  terms  will  be  discussed  in  the  following  paragraphs. 

1)  The  heat  transfer  into  the  parachute  by  convection.  Q^,  is  cal¬ 
culated  from  the  relation 


hS, 


(T  -  T  ) 
aw  w 


(9) 


where 

«  total  surface  area  of  decelerator  exposed  to  convective 
heating,  (ft^) 

o 

h  «  convective  heat  transfer  coefficient,  (BTU/f t*’-sec-°R) 
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T  =  adiabatic  wall  temperature,  (°R) 
aw 

T  =  temperature  of  the  surface  of  the  dec.elerator ,  (°R) 

The  flat  plate  with  constant  wall  temperature  serves  as  the  standard 
condition  to  which  simple  formulae  apply.  For  this  case  the  differential 
equations  which  describe  the  flow  and  energy  conditions  are  well-estab¬ 
lished  for  a  laminar  boundary  layer  at  normal  conditions  of  pressure 
(continuum  flow  regime)  and  temperature  (ideal  gas  behavior).  The  prin¬ 
cipal  difficulties  in  solving  these  equations  stem  from  the  large  varia¬ 
tions  with  temperature  and  pressure  of  such  thermal  properties  as  vis¬ 
cosity,  thermal  conductivity,  density,  and  specific  heats.  Experimental 
imperfections  generally  lead  the  designer  to  conclude  that  in  the  laminar 
range,  the  most  reliable  information  is  obtained  from  analytical  calcula¬ 
tions.  If  the  geometry  is  complex  and/or  transition  to  turbulent  flow 
takes  place  on  the  body,  then  empirical  methods  must  be  resorted  to. 

The  proposed  methods  for  the  calculation  of  h  are  summarized  in  Section  II. 

2)  The  energy  gained  due  to  incident  solar  radiation  may  be  deter¬ 
mined  from  the  relation 


las 


so 


rad 


(10) 


whe  r  e 

1  -  solar  radiation  intensity;  a  reasonable  estimate  of  its  value 
is  425  BTU/ftZ-hr 

a  “  absorptivity  for  solar  radiation  of  the  parachute  material;  a 
typical  value  is  0.2 

2 

S  *°  ProJecve^  surface  area  of  body,  (ft  ) 
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3)  The  energy  loss  from  a  parachute,  by  radiation  to  space  is  given  by 


sp 


aS 


rad 


F 

p-sp 


T*  ) 
space 


(ID 


-8  2  4 

a  *  Boltzmann  constant  equal  to  0.173  x  10  BTU/ft  -hr-°R 
Fp  Sp  *  radiation  shape  factor  between  parachute  and  space 

e  **  emissivity  of  the  parachute  material  for  long  wavelength 
radiation;  a  typical  value  is  0.8 
T  ■  temperature  of  the  walls,  (°R) 

T^pace  “  temperature  of  space,  (7°R) 

4)  The  energy  loss  front  the  parachute  to  the  earth  is  given  by 
the  relation 


where 


"Vrad  £  Vearth  '  Te,rth> 


(12) 


Fp  eart^  *  radiation  shape  factor  between  the  parachute  and  the  earth 
e  **  emissivity*  has  the  same  significance  and  practically 


the  same  value  as  in  Equation  (11) 


T  arth  "*  average  earth  temperature;  assumed  to  be  520“R 


5)  The  rate  at  which  energy  is  stored  by  the  parachute  is  given 


by  the  relation 


'at 


p  V  c 
m  m  p 


dT 

dT 


(13) 
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where 


=  average  density  of  parachute  material  (lbm/ft  ) 

3 

V  =  total  volume  of  parachute  material  (ft  ) 
m 

c  =  heat  capacity  of  parachute  material  (BTU/lbm-°R) 

Pm 

d” 

■—  =  time  rate  of  change  of  the  average  temperature  f  the  para¬ 
chute  mass,  (°R/hr) 


Total  Enerf  valance 

An  energy  balance  may  now  be  written  based  on  the  abov^  equations 
to  yield 

hS  (T  -  T  )  +  IaS„  -  oS,t(T4  -  T4  )  + 

1  aw  w  2  3  w  space 


aF.S,  c(TA  -  T4  } 
h  4  w  earth 


„  dT 
p  V  c  -r~ 
m  m  p  at 
r  m 


(14) 


In  a  majority  c*f  cases  considered  the  terms  for  the  energy  flow  into  the 
parachute  by  solar  radiation  and  for  energy  loss  due  to  radiation  to 
the  earth  are  generally  negligible  compared  to  those  of  radiation  to 
space.  In  addition  the  temperature  of  space  (Tg  ,c  ■'  7°R)  may  also 
be  neglected  compared  to  equilibrium  temperatures  of  practical  interest. 
Finally,  if  no  heat  is  absorbed  by  the  parachute  the  convective  flux 
must  be  equal  to  the  radiant  flux.  This  condition  permits  a  simple 
determination  of  the  equilibrium  surface  temperature  T  .  Consequently 


hS,  (T  -  T  )  -  oS  f.T 
j  aw  e  3  e 


23 


(15) 


or  rewritten  as 


T  -  T 
aw  e 


(if) 


Equation  (16)  may  be  solved  for  the  equilibrium  temperature  in  the  above 
form  by  trial  and  error >  or  the  following  system  may  be  used  which  is 
an  explicit  solution  for  Tp  from  Equation  (16): 


T  -  ~ 
e  2 


— y'L  +  /L  -  2(1.  -  m/rL) 


where 


li®!.  +  >45  + 

j_  2  V  4  27 

tft)  •  -frfe) 


1/3 


aw 


(17) 


(18) 


(19) 


Opening  Period 

Thv  extremely  short  time  to  maximum  temperature  of  sotre  parachutes 
should  be  noted  since  it  becomes  of  the  same  order  of  magnitude  as  the 
opening  time  (which  falls  in  the  range  of  1/20  to  1/2  second,  depending 
on  thu  parachute  design).  It  is  estimated  in  Reference  10  that  the  time 
to  maximur  temperature  should  be  at  least  ten  times  greater  than  the 
parachute  opening  time  before  the  opening  effect  can  be  neglected. 

The  greatest  source  of  uncertainty  in  Equation  x16)  lies  in  the 
evaluation  of  h.  It  has  been  demonstrated  that  the  dimensionless  heat 
transfer  coefficient  (the  Mussel t  number*  for  the  constant  property,  in¬ 
compressible  flow  case  may  be  utilized  to  ca’culate  heat  transfer  for 
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Mill 


the  variable  property,  supersonic  case  if  the  various  physical  properties 
which  appear  are  evaluated  at  a  reference  temperature  T  .  According  to 
Reference  11, 

T  -  1'  +  0. 5(T  -  T,)  +  0.2><T  ~  T„)  (20) 

r  6  w  6  aw  6 

The  recovery  temperature  is  dependent  on  the  trajectory  and  will  vary 
with  altitude.  At  a  given  altitude  the  recovery  temperature  is  that 
temperature  which  a  parachute  would  ultimately  assume  under  the  sole 
influence  of  steady  state  convective  heat  transfer,  i.e.,  in  the  absence 
of  other  modes  of  heat  transfer.  The  recovery  temperature  is  determined 
by  the  extent  the  kinetic  energy  of  the  air  stream  is  converted  into 
thermal  energy.  For  blunt  shapes,  the  conversion  is  very  nearly  complete 
such  that  the  recovery  temperature  is  closely  equal  to  the  stagnation 
temperature. 

According  to  Equation  (20),  the  calculation  of  the  reference  tem¬ 
perature  (T  )  requires  advance  knowledge  of  the  unknown  wall  temperature. 

For  blunt  bodies,  normally  Tr  -  T  >>  T  and  T  =  i/2  (T„)  such  that 

o  aw  w  r  0 

the  reference  temperature  T  is  insensitive  to  T  for  the  early  portion 

r  w  J 

of  the  trajectory. 
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SECTION  III.  CONVECTIVE  HEAT  TRANSFER  COEFFICIENT 


In  the  selection  of  the  film  heat  transfer  coefficient  it  is  necessary 
to  determine  the  characteristics  of  the  viscous  flow  field  adjacent  to 
the  surface  elements  of  the  decelerator.  There  are  no  exact  carcula- 
ticns  of  "supercritical  woven  me  h  flew"  so  that  2  rigorous  development 
is  not  possxble.  There  have  been  two  analytical  models  proposed  for  the 
prediction  of  the  average  convective  heat  transfer  coefficient  for  mesh 
flows  of  intermediate  porosity  (single  bow  shock-individual  boundary 
layers  on  each  filament).  They  are  as  follows: 

1;  two-dimensional  laminar  stagnation  flow  coupled  with  stagnation 
pci  t  velocity  gradients  for  low  speed,  potential  flow  over  an 
infinite  cylinder 

and 

2)  axisymmetric  turbulent  boundary  layer  flow  as  applied  to  the 
throat  region  of  a  supersonic  nozzle 

Rose,  Prof  stein,  and  Adams  (12)  have  reported  that  the  critical 
Reynolds  n umbel  for  the  transition  from  3  laminar  to  a  turbulent  boundary 
I  aver  is  approximately  4  a.  10^,  based  on  conditions  at  the  outer  edge 
of  the  boundary  layer. +  With  fiber  diameters  on  the  order  of  0.004  -  0.050 
inch  xgure  4  reveals  that  a  laminar  boundary  layer  will  exist  over  the 

A  somewhat  mere  exact  and  complex  criterion  for  transition  is  that 
of  an  empirically  determined  critical  Reynolds  number  brsed  not  on  the 
wetted  length  but  <\n  the  momentum  thickness.  The  critical  momentum  thick¬ 
ness  Reynolds  number  is  approximately  250  (13). 
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entire  front  face  of  a  typical  mesh  element  for  all  Mach  numbers  and 
altitudes  under  consideration  for  woven  mesh  parachutes. 

The  experimental  results  of  Scott  and  Ruiz-Urbieta  (14;  Figures  5-8) 
produced  an  engineering  conclusion  that  there  are  no  large  variations  of 
local  heat  flux  around  the  periphery  of  a  woven  mesh  surface  and  that 
the  conventional  laminar  stagnation  point  prediction  is  a  useful  first 
approximation.  In  this  paper,  local  measurements  were  made  on  an  en¬ 
larged  woven  mesh  in  which  the  individual  fibers  were  one  inch  in  diameter. 
Measurements  of  the  local  pressures  (Figure  6) ,  recovery  temperatures 
(Figure  7),  and  convective  heat  transfer  coefficients  (Figure  8)  were 
made  on  the  forward  and  rearward  surfaces  of  the  mesh  at  locations  in¬ 
dicated  in  Figure  5.  Pressure  ratios  across  the  mesh  were  varied  from 
subcritical,  through  the  critical,  and  into  the  supercritical  regime. 
Comparisons  of  the  experimental  data  with  the  two  proposed  methods  of 
calculating  average  neat  transfer  coefficients  are  readily  observed. 

It  is  of  importance  that  no  large  variations  of  h  appear — a  result  which 
justifies  the  use  of  a  single  heat  transfer  coefficient  applied  to  tne 
entire  surface.  It  is  of  interest  here  that  the  experiments  revealed 
that  regions  of  separated  fl^~  covered  approximately  70  percent  of 
the  surface  area  exposed  to  convective  heating.  Secondly,  comparisons 
of  the  data  with  both  predictive  models  mentioned  above  revealed  that 
while  the  two  methods  were  generally  conservative,  some  of  the  data  taken 
at  locations  between  the  stagnation  line  and  the  sonic  throat,  do  fall 
above  the  stagnation  point  predictions.  Since  the  experimental  model  was 


Figure  6.  Pressure  Distributions  on  Individual  Mesh  Elements 
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Cross-Section 


Figure  6.  Heat  Transfer  Coefficients  on  Mesh  Elements  and  Comparisons  with  Analyses 


hollow  and  of  thin-walled  construction,  internal  heat  conduction  could 
not  damp  out  these  local  variatir  is — as  would  a  solid  model. 

It  is  concluded  that  the  laminar  stagnation  point  approximation 
should  be  used  in  the  predictive  scheme  with  the  turbulent  calculation 
also  being  carried  out  to  indicate  an  upper  (extreme)  limit  on  the  con¬ 
vective  heating  process. 

Hie  discussion  is  also  restricted  to  heat  transfer  under  steady 
state  conditions.  It  has  been  shown  that  for  rates  of  change  in  velocity 
as  they  are  usually  encountered  in  aircraft  the  convective  heat  transfer 
may  be  regarded  as  quasi-steady.  This  means  the  heat  transfer  at  each 
instant  may  be  calculated  with  formulae  valid  for  the  steady  state  into 
which  property  values  are  introduced  as  they  are  encountered  at  this 
specific  instant. 


Laminar  Stagnation  Point  Heat  Transfer  Coefficient 

In  the  neighborhood  of  a  stagnation  point  where  the  velocity  dis¬ 


tribution  is  given  by  Ux 

6 


Cx,  the  heat  transfer  coefficient  is  given  by 


(21) 


if  the  heat  generation  due  to  compressive  work  and  viscous  dissipation 

is  neglected.  The  values  of  A  given  by  Goldstein  (15)  depend  upon  the 

Prandtl  number,  Pr.  A  few  values  of  A  are  listed  in  Table  I. 

Table  I.  The  constant  A  for  a  circular  cylinder 
heat  transfer  coefficient 


0 . 7 _ 0. 8  0.9 _ JyO _ y 

0.495  0. 521  0.546  0.570  0.592 
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1 

For  a  circular  cylinder  (x)  =  211^8^ (x/tjID)  and,  therefore,  C  = 

(dU^/dx)^  ^  ~  411^/D  in  the  neighborhood  of  the  stagnation  point.  Con¬ 
sequently  the  Nusselt  number  based  on  the  diameter  of  the  cylinder  becomes 


Nu  » 


(22) 


Combined  Analysis  for  the  Heat  Transfer  Coefficient 
Referring  to  Figure  2  and  Equation  (22),  the  principal  problem  is  to 
calculate  the  equivalent  approach  flow  velocity  for  the  isolated,  in¬ 
finite  cylinder. 

By  continuity 


P3V3A3  "  P*a*A* 


(23) 


The  geometric  porosity  is  defined  as 


P  -  A*/ A 


total 


(24) 


and  Equation  (23)  is  written  as 


P3V3  _  A* 
p*a* 


total 


total 

A. 


total 


The  geometric  factor  Atota]/A3  determined  from  the  parachute  cross 
section.  We  consider  the  flow  entering  the  area  A^  to  originate  in  a  fic¬ 
titious  reservoir  where  the  temperature  is  T^,  pressure  p'q.  and  density  p'q* 
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Since 


p*a*  -  pA  ^RT*  «  -®T  /jP  o'  ft  ( 25 ; 

p0  p0  V  0  u  o 

We  shall  assume,  for  the  purpose  of  the  heat  transfer  calculation  only, 
that  =  Pq"  This  will  prove  to  be  conservative  as  far  as  surface 
temperatures  are  concerned. 

Equations  (24)  and  f?5)  are  combined  to  yield 


/E  _ 

[p6  h 0  *3  P  Atotal 


Tha  velocity  v^  is  taken  to  represent  the  approach  flow  velocity 
required  in  the  laminar  stagnation  point  analysis.  Equation  (2o)  is 


rewritten 


Uoo  "  v3  * 


total  Up*  A*  r 


Referring  to  Figure  9,  the  porosity  P  is  given  by 


&  2 
*  total  { u  +  D) 


where  L  is  the  distance  between  the  i  i lament,,  and  I)  is  the  diameter  of 
any  filament .  Therefore  the  expression  for  He  approach  velocity  becomes 


(L  +  0) 


>)Hk  M  * ' 


Typical  Ro j :  Nozzle  Geometry 


5 

i 
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For  a  Frandtl  number  of  0.71,  the  Goldstein  result  (Equation  [22 ])  for 
the  heat  transfer  coefficient  to  an  infinite  cylinder  is 


h 


0.496 
Pr  ' 


(30) 


Now  with  pi  ■  p!  8  i  ,  *  d  .  8  T  ,  p  «  n'  because  of  the 
*  0  o  0  ref  ref  ref  *  *ref  *T> 

condition  of  a  constant  static,  pressure  across  the  stagnation  point  boun¬ 
dary  layer. 


If  we  make  the  relatively  common  assumption  that  p  =  CT,  that  is 


P0 
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the  expression  for  the  heat  transfer  coefficient  becomes 


leading  sonic  throat 

constant 


geo-  wire  shape  trejet^ory 

metric  size  factor  variables 

weave 


(33) 


The  heat  transfer  results  for  a  variety  of  geometries  arc  an  If ied  through 
the  use  of  the  sonic  Reynolds  number  Re*.  With  the  aid  of  equation  (2b) f 
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r 


Re*  «  p*a* 


l:* 


p* 
_  » 


”i/x 
T*  VR 


D 


■*o 


(34) 


and  the  final  relation  between  h  and  Re*  becomes 


Nu*  =  ~  -  2(0.496) 


(-Hfc  ®  ®  (35) 


Turbulent  Analysis 

Sibulkin  (16)  developed  an  approximate  solution  for  the  heat  t  ansfer 
to  a  constant-property,  turbulent  boundary  layer  with  arbitrary  freestream 
velocity  and  surface  temperature  distributions.  In  order  to  estimate 
rapidly  the  maximum  rate  of  heat  transfer  which  occurs  near  the  throat 
of  a  supersonic  nozzle,  a  throat  approximation  was  developed  which  involves 
only  the  nozzle  supply  conditions,  throat  opening,  and  the  throat  radius 
of  curvature.  The  Sibulkin  prediction  leads  to  the  following  relation 
for  h: 


where 

2 

h  =  heat  transfer  coefficient,  (BTU/sec~f t^-°R) 

C  =  axially  symmetric  constant,  (C  *  0.0027  BTL'/ lb-sec^^-f t  ^ "*) 
r*  *  nozzle  throat  radius,  (ft) 

L*  =  radius  of  curvature  tt  nozzle  throat,  (ft) 
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Pq  «*  nozzle  supply  pressure,  (lb^/ft  ) 
u*  *  dynamic  viscosity  at  T,  ( slugs/f t-sec) 

p*  *■  sonic,  density  at  T,  (slugs/ft^)  =  (0.5283  p^)/ (1716) (0.833  Tq)-Tq/T 

Tq  »  nozzle  supply  temperature,  (°R) 

T*  **  sonic  temperature  =  0.833  Tq,  (°R) 

I  "  (T*  +  T  )/2,  reference  temperature,  (°R) 

I  w 

The  velocity  gradient  at  the  nozzle  throat  of  an  axisymmetrical-f low 
nozzle  is 


The  assumed  circular  nozzle  throat  of  radius  r*  is  frequently  only  an 

approximation  to  the  actual  rectangular  geometry.  That  is,  the  orifice 

2 

shown  as  the  shaded  area  in  Figure  9,  r.r*  ,  is  assumed  to  be  circular 

2 

with  the  same  geometric  area  as  the  actual  rectangular  orifice,  L  .  The 
corresponding  stagnation  point  velocity  gradient  on  an  infinite,  isolated 
cylinder,  reduces  to 
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SECTION  IV.  INTERNAL  ENERGY  TRANSFER 


The  heat  gained  by  a  body  is  determined  by  the  resistance  to  the 

flow  of  heat  at  the  surface  as  well  as  the  internal  resistance.  Consider 

a  wall  losing  heat  by  convection  to  an  ambient  gas.  At  the  wall,  the 

convective  velocity  is  zero  such  that  in  the  absence  of  radiation,  the 

energy  transport  is  by  conduction  only.  The  Fourier  conduction  law 

states  that  a  rate  of  heat  flow  dQ  in  the  direction  of  the  normal  n  to 

an  area  element  dA  is  dQ  =  -  kdA  3T/3n.  The  minus  sign  indicates  that 

heat  flows  in  the  direction  of  decreasing  temperature.  According  to  the 

Langmuir  film  concept  the  temperature  gradient  in  the  gas  is  confined  to 

a  narrow  layer  of  thickness  6.  We  simplify  the  convective  heat  exchange 

process  in  the  gas  by  assuming  the  entire  thermal  resistance  to  the  heat 

exchange  between  the  wall  and  the  gas  is  concentrated  in  a  static  film 

which  adheres  to  the  wall,  whereas  outside  this  film  all  temperature 

differences  vanish  as  a  result  }f  the  mixing  motions  of  the  gas.  The 

flow  of  heat  to  the  wall  is  given  by  =  Qc^/6^)  A  (T^  -  T^) .  Owing  to 

the  difficulty  in  interpreting  the  film  thickness  6  ,  it  is  common  en- 

g 

gineering  practice  to  combine  the  quotient  k  /6  into  a  single  quantity 

g  g 

a  isignated  h,  the  film  heat  transfer  coefficient,  and  express  the  boun¬ 
dary  heat  exchange  as  Q  »  h  A  (T  -  T  )  where  A  is  the  boundary  area  in- 

g  w 

volved.  The  quantity  hA  is  referred  to  as  the  surface  conductance  while 

the  value  1/hA  =  R  's  called  the  surface  thermal  resistance  of  the  film 
3 

heat  transfer  process.  For  combined  excharge  processes  the  sum  of  the 
individual  resistances  is  the  thermal  resistance  of  the  overall  heat 
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transfer  process.  The  dimensionless  convective  heat  transfer  coeffi¬ 


cient  is  called  the  Nussclt  number  Nu  *  h6  /k  and  assumes  values  both 

conv  g  g 

greater  than  and  less  th :n  the  value  of  unity  predicted  by  the  oversim¬ 
plified  static  film  concept. 

The  heat  entering  the  body  is  also  governed  by  the  Fourier  conduction 

law.  Whenever  the  internal  conduction  is  of  importance,  a  dimensionless 

conduction  Nusselt  number,  given  by  hu  ,  =  h6  /k  **  Bi  (usually  referred 

°  cond  w  w 

to  as  the  Biot  number),  enters  as  a  parameter.  The  Biot  number  can  be 
thought  of  as  the  ratio  of  internal  to  surface  resistance  since 


Bi 


hi 

k 


(38) 


It  is  of  primary  importance  to  estimate  the  value  of  the  Bio^  number  in 
a  given  problem  in  order  to  justify  any  further  assumptions  used  to 
simplify  the  analysis. 

In  any  problem  of  transient  heating  (cr  cooling)  the  heat  transfer 
process  is  influenced  by  both  the  internal  resistance  R^  ,  and  the  surface 
resistance  R  .  The  limiting  cases  are  represented  by  a  zero  internal  re- 
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Case  A  —  *  0,  Negligible  Internal  Resistance 

When  the  body  is  constructed  of  material  possessing  a  large  thermal 
conductivity,  its  internal  thermal  resistance  may  be  ignored  so  that  the 
heat  transfer  process  is  controlled  only  by  the  surface  resistance.  No 
temperature  differences  can  exist  in  the  body  under  these  conditions- -the 
body  is  isothermal.  Consider  a  body  having  an  initial  temperature 
at  t  =  0  that  is  suddenly  exposed  to  a  fluid  at  temperature  T^.  The  rate 
of  energy  transfer  from  the  surroundings  for  t  >  0  is  given  by 


Q  -  p  V  c  ~  =  hA  (T,  -  T) 
x  m  m  p  dt  f 

m 


(39? 


where  c  ,  p  ,  V  ,  and  A  are  the  specific  heat  capacity,  density,  object 

m  _  _  __ 

volume,  and  object  surface  area,  while  T  =  T(x)  ■  T  (t)  is  the  uniform 

w 

temperature  at  any  time  x  and  h  and  T^.  are  the  heat  transfer  coefficient 
and  ambient  fluid  temperature.  We  shall  eliminate  the  bars  for  con¬ 
venience.  If  h  and  are  Independent  of  x.  Equation  (39)  integrates 
to  yield 


/hA - \  t 

\  o  c  V  ) 

'  m  p  m  / 

*  m 


T,  -  T 


(AO) 


f 


and  the  uniform  temperature  of  the  object  rises  exponentially  with  time. 
If  we  introduce  the  Biot  number  ii  i  and  Fourier  number  Fo  defined  by  the 
relat ions 


JT 


k  i 
w 


w 


iM 


Fo 


(41) 


r  or 


Convect-  j  Environment  t  at  r  *  R-,  High 


ure  13.  (Continued 
(Reference  .20) 


Equation  (40)  becomes 


T  -  T, 


Ti  *  Tf 


oiFo 


(42) 


The  product  BiFo  is  independent  of  the  object  thermal  conductivity  k  . 

W 


For  infinite  cylin<‘«rs  the  geometric  factor  (A6  _/V)  =  2  and  Equation  (42) 
is  rewritten 


2hr 


T  -  Tf 

T  -  T 
i  f 


pc  6 

m  p  w 
m 


(43) 


If  the  temperature  of  the  external  fluid  varies  linearly  so  that  T 


f 


_I  C  (t),  the  solution  becomes 
*i 


T  -  T  +  C 


p  c  o 

m  p  w 
m 

2h 


2hr 


pc  6 

m  p  w 
m 


Tt  -  Tf  +  C 


pc  P 

ran  w 
'  m _ 

2h 


(44) 


-C  i 


If  the  fluid  temperature  diminishes  exponentially  so  tuat  T .  »  i\ 
the  solution  bvLuttir»  t.  i  /  ) 


111! 


(  T 


(I, 


T  j  -  T  (,  -  1) 
___  ______ 

i  l 


m  p  w 
m 


(  4  b ) 
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where 


P 


1  - 


-  1 


c  c  6  C 
m  p  w 
m 


2h 


(46) 


If  the  temperature  of  the  external  medium  varies  periodically  (one-half 
cycle  is  an  approximation  to  the  atmospheric  entry  case)  such  that  * 

+  C  cos  (air),  where  T  is  the  mean  temperature,  of  the  fluid,  constant  with 
time,  and  S  and  w  are  the  amplitude  and  frequency  of  oscillation  of  the 
temperature  difference  the  solution  is 


2’n 


pc  5 
m  p  w 

T  -  I  -  A'  m 

- — - - —  =  e 

T.  -  T  -  A 
l 


(47) 


where 


A’ 


cos  wt  +  a)  sin  u>x 

- 


P 


2h 


pc  6 

m  p  w 
m 


2 
ij)  ’ 


(48) 


w.-isi'  li  ~  -  K  *  '),  Ncgl  igib  1 e  Surface  Resistance 

A  class  of  problems  exists  in  which  the  surface  thermal  resistance 
is  negligible  and  the  internal  thermal  resistance  dnminat.ee  the  problem. 
If  the  ambient  temperature  of  the  surrounding  fluid  is  T  (i)  we  take  the 
temperature  of  the  surface  of  the  object  to  be  T  ( i )  for  i  U. 
a.)  T.  -  constant  T.  =  T;(r) 

Consider  an  infinitely  long  circular  cylinder  with  the  initial 
temperature  distribution  T,(r,0},  The  temperature  of  the  surface  T{R,i) 
is  suddenly  changed  to  T  --  const  for  all  r  >  0.  The  solution  roust  sat¬ 
isfy  die  system 


92T  1  _3T 

^  2  l  O  L 

c  5T 


(pc  ) 
_ 

K. 


3T 

3t 


_i  _di 

<*’  O  t 


(49) 


T  = 

1\  (r) 

a  t  !  -•  0  ; 

0  <  r 

T  = 

0 

at  r  =  R; 

0  •>  0 

(50) 


For  the  special  case  of  a  uniform  temperature  distribution  T  (r) 


V 


the  solution  of  the  above  system  is 


T  -  Tf 
Ti  -  Tf 


«  2  V 


M*Fo  ,  /M  r) 
n  J0  (  Mn  K  ) 


i  M 

n=l  n 


J,(M  ) 
1  n 


(51) 


where  J,.  and  J  are  the  zeroeth  and  first -order  Bessel  functions  and  M 
u  l  n 


represents  the  roots  of  the  zero-order  Bessel  function  J„(M  )  =  0. 

(J  n 


i  i  rs  t 

‘  ive  roots  ol  . 

i  (M  )  are 
m  n 

>;,iv<.*u  i 

U  tae  : 

lull  owing  table 

(  i  8 .)  : 

Tab  1  e 

IJ.  Hoot; 

:•  M  O  J 
n 

.1  (H  ) 

m  u 

■'  0 . 

m 

Mi 

M2 

M3 

M, 

a 

M5 

0 

2.4048 

r  .  5201 

8.6 

5  37 

I  J  .  79  1 .5 

14.9 309 

1 

3.8317 

7 . 0156 

10. 1 735 

1 3.3237 

16.4706 

The  series  lias  been  computed  for  the  axial  temperature  history  at.  the 


centerline  r  *  0  according  to  the  following  expansion: 
T (0 ,  t  )  ~  Tf 


A  few  values  of  C(Fo)  are  given  in  the  following  table: 

Table  III.  Values  of  C<  Fo)  in  Equation  (52) 


-  M^Fo 
n 

Of 

-  M,  Fo  -  xM^Fo 

V  e _ 

42 _  ^  _ _ _ __  ^ 

L,  H  J,(M  ) 
n=l  n  J  n 

*“  z 

M.J.frJ.)  M0 J  .  (M,, ) 

i  l  tL  *  A. 

Fu ,/ 


(52) 


Fo 

C  (  Fo ) 

0.020 

0.99999 

0.040 

.99627 

0.060 

.97054 

77? 

0.100 

.848 36 

0.120 

. 77293 

0.  140 

. 69798 

0.  160 

. 62692 

Fo  C  (  Fo  ) 

0.200  .50149 

0.240  . 39912 

0.280  .2.1704 

n  ,'OH  75)67 

0.360  . 19973 

0.400  .  1585 

0.500  .0887 

1.000  .0049 

i  .  500 


5i 


0.180 


,  56  I 2n 


.000  3 


constant 
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Cast-  C  --  KJ n'i  If  Internal  and  Sur  face  Resistance 


The  most  general  heat  transfer  problem  involves  heating  or  cooling 
of  an  object  in  which  both  the  internal  and  surface  resistance  must  he 
considered.  Let  the  initial  temperature  throughou*  -finitely  long 

cylinder  be  ’’h  =  constant.  At  the  time  t  =  0  the  cylinder  is  exposed  to 
convective  heating  from  an  ambient  fluid  at  Tj.  =  constant .  In  addition 


the  heat  transfer  coefficient  is  taken  to  be  constant. 


.he  solution  .is 


T  -  T  ”  1 

_ 1  ..  a  V 

T .  -  T 1  L  ■' 
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n=l  Mn  fjf  (M  )  +  J^M  )~| 
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,(%  t) 


(56) 


in  which  M  are  the  roots  of  tne  eigenfunction  Bessel  equation 

J,  (M  ) 

M  T~TZTT  “  Hi 
n  J0(Mn) 


(57) 


Note  thstt  the  temperature-position  solution  is  a  function  of  two  parair;— 
eters  (Fo  =  ax/R^,  Bi  =  hR/k) ,  rather  than  one  as  in  tiie  case  of  negli¬ 
gible  surface  resistance.  Tne  first  three  roots  of  Equation  (57)  are 
given  in  the  following  table: 


Bi 

M2 

M 

“3 

0. 20 

0.6J  70 

3.8835 

7 . 0440 

0„  40 

0.8516 

3.9344 

7.0723 

O 

o 

rH 

1.2558 

4.0795 

7.1558 

5.00 

1.6898 

4.7131 

7.6177 

10.00 

2,1795 

5.0332 

7.9569 

100.00 

2 , 3809 

5.4652 

8.5678 

For  values  of  Fo  greater  than  0.2,  all  terms  In  Equation  (56;  after  the 
first  become  negligible.  Therefore,  for  r  equal  to  zero 

-  M?Fo 

KO)  -  Tf  2J1(M1)e 

Ti  "  Tf  M1[Jy(M1)  +  jjCMj)] 

With  this  approximation  in  mind,  Heisler  has  presented  time-temperature 

charts  of  the  solution  of  Equation  (56)  for  the  range  of  Biot  uuw>.ut 

from  Bi  *  0.01  to  100.  The  Heisler  charts  employ  Equation  (58)  for  the 

central  temoeratn h*st<">rv  =  i.01  it  is  clear  that  for  any 

u 

other  position  0  < r  <  R,  (T(r)  -  Tf)/(Ti  -  Tf)  -  (T(0)  -  Tf )/(T±  -  T  ) 

•  f  J .  A  few  useful  values  of  the  zeroeth  and  first-order  Bessel 

0  \  r  K/ 

functions  of  the  first  kind,  J_(M)  and  J. (M  ),  are 

U  In 
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Tab  1 1-  V. 

M 

n 

Zu.roeth  and  First -Order  B< 

d,.  (M  ) 

0  n 

.'skim  Fount  ions 

J,  (M  ) 

1  n 

0.0 

1.0000 

0.0000 

0.2 

0.9900 

0.0995 

0.4 

0.9604 

0.1960 

1.0 

0. 7652 

0.4400 

1 . 5 

0.5118 

0.5579 

2.0 

0.2239 

0.5767 

2.5 

-0.0484 

0.4971 

It  is  that  if  Bi  <  0.1,  <  0.5  and  J^(0.5  r/R)  >  0.95.  It  is 

see"  that  the  temperature  distribution  through  the  cylinder  is  nearly 
uniform  for  Bi  <  0.1.  Therefore,  for  Bi  <  0.1  the  combined  heating  process 
may  be  considered  to  be  controlled  by  the  surface  resistance  and  the  body 
uniform  temperature  is  governed  by  the  Newtonian  equations. 

Fox  very  large  values  of  Biot  number  or  for  small  values  of  the  Fourier 
number  the  series  solution  given  by  Equation  (56)  converges  slowly.  For 
this  reason  Heisler  developed  "short-time"  charts  for  1/Bi  *  0.G5  to  00 
and  Fo  =  0  to  0.2.  For  Bi  =  10,  (T(R)  -  Tj.)/(1\  -  T^.)  reaches  the  value 
0.5  when  Fo  =  0.01.  In  this  case,  the  internal  resistance  prevents  the 
passage  of  heat  into  the  interior  so  that  tne  surface  layers  "'ust  absorb 
the  energy  input.  The  small  mass  available  leads  to  rapid  temperature 
increases  near  the  surface  such  that  the  negligible  surface  resistance  case 
is  approached  in  which  the  surface  temperature  instantaneously  assumes  the 
fluid  temperature. 
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/he  asymptotic  expansions  for  small  values  of  time  are  cumbersome 
and  an-  valid  ciiiy  j  f  r/R  is  not  small.  Heisler  computed  the  short -tim 
curves  lor  the  cylinder  surf  acr  while  the  other  charts  were  obtained  by 
"electrical  analogy"  on  a  "heat  and  mass  flow  analyzer."  Charts  for 
r/R  “  0,  U.5f  and  1.0  are  presented  as  Figures  14  and  15. 

Case  D  Approximate  Methods 

Of  course  the  only  method  for  solution  ot"  the  general  problem 
(0  <  Bi  <  with  variable  boundary  conditions  [h(i),  T  (x) ]  involves 

3W 

numerical  integration  of  tne  heat  conduction  equation.  Short  of  this, 
is  often  useful  to  obtain  rapid  esffkaates  of  the  surface  temperature,  1 
using  the  results  of  Sections  IV.  Case  A  and  IV,  Case  B  to  predict  the 
results  of  Section  IV,  Case  C.  Consider  the  two  sketches  below: 


The  problem  is  to  predict  the  intermediate  Bi  behavior  from  the  limiting 
cases  as  is  illustrated  by  the  dashed  lines  in  the  sketches  above.  Asymp¬ 
totic  expansions  in  the  two  ranges  Bi  >>  0,  Bi  -  0  are  reasonable.  A 
somewhat  more  useful;  though  cruders  approximation  is  the  following: 


T  -  T 
Bi  Bi-0 

T  -  T 
Bi*®  Bi-0 


tanh 


where  C  is  an  arbitrary  constant.  The  hyperbolic  tangent  was  chosen  be¬ 
cause  it  has  the  proper  asymptotic  behavior  (tanh  0  -  0,  tanh  <*  -  1)  . 
Below  is  a  short  table  of  the  hyperbolic  tangent. 

Table  VI.  Table  of  Tanh  x 


X 

tanh  x 

X 

tanh  x 

0 

.0000 

1.2 

.8337 

0.2 

.1974 

1.4 

.8854 

0.4 

.3800 

2.0 

.9640 

0.6 

.5370 

3.0 

.9951 

0.8 

.6644 

4.0 

.9993 

1.0 

.7617 

5.0 

.9999 

(If  x  >  5,  tanh  x  *»  1.0000  to  four  decimal  places.) 
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A  sample  calculation  may  be  carried  out  using  the  Heisier  charts  (Figures 


in  and  17).  For  Fo  =  0.2,  Figure  ll>  reveals  the  following  data: 


El. 

T(R)  -  T 

tanh  (fbT/2) 

T  -  1  . 
g  i 

1 . 00 

1 . 00 

10 

0.92 

0.92 

0 . 84 

0.81 

1.00 

0.43 

0.46 

.1-0 

0.26 

0.33 

.10 

0.06 

0.13 

0 

0 

0 

The  approximation  is  examined  for  C  =  2  above  an -i  is  reasonable  in  the 
range  0.5  <  3i  <  10.  This  is  the  useful  range  of  most  parachute  materials. 
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INFINITE 

CYLINDERS 


I 


8| 


Figure  17.  Heisler  Position  Correction  Factors  for  Dimensionless 
Temperature  Ratios  for  Infinitely  Long  Cylinder. 

To  be  Used  with  Figure  16  (Reference  21) 
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APPENDIX  I. 


APPLICATION  OF  THE  PROPOSED 
TECHNIQUES  TO  A  PRACTICAL  EXAMPLE 


The  analytical  tools  developed  in  the  text  were  applied  to  a  para¬ 
chute  which  was  launched  by  the  Goodyear  Corp.  (22).  The  heat  transfer 
calculations  are  based  on  measured  values  of  velocity  versus  altitude 
(Figure  35  in  Ref.  22)  and  calculated  values  of  the  Mach  number,  total 
pressure,  and  total  temperature.  The  important  aerothermodynamic  vari¬ 
ables  are  tabulated  in  Tables  VII  and  VUL  .  The  data  required  for  the 
calculation  of  the  convective  heat  transfer  coefficients  (Equations  [33] 
and  [35])  are  presented  as  the  last  two  columns  on  the  right  in  Table  VIII 
These  date  are  plotted  in  Figures  IB  through  20.  It  is  observed  that 
the  parachute  was  deployed  during  the  ascent  phase  of  the  trajectory. 

Experimentally  (23)  it  has  been  found  that  the  actual  pressure  behind 
the  normal  shock  is  less  than  p^  calculated  using  normal  shock  theory  in 

the  absence  of  any  forebody  wake.  The  ratio  p^  /p^  is 

actual/  calculated 

given  in  Reference  23  for  two  Mach  numbers,  i.e.. 


M 

M 


1 

1 


3 

4 


*  0.61 


-  0.47 


Since  the  present  data  fall  in  the  range  4  <  <  5.5  extrapolation  was 

necessary.  The  following  asymptotic  behavior  seems  plausible; 


M 


1 


1 


1 


a 


It  is  obvious  to  try  the  forms 


I  (M£  -  1) 


or  e 


B(MX  -  1) 


The  best  fit  gave 


M1  -  1 


The  general  heat  balance,  Equation  (14),  was  applied  to  one  mesh,  silent" 

ly  assuming  that  by  rotation  and  fluctuation  of  the  parachute  the  local 

variations  in  hearing  due  to  location  on  the  parachute  are  leveled  out. 

From  the  given  porosity  P  ®  0.13  and  the  cross-section  of  one  wire 

(see  sketch  below)  one  can  calculate  the  mesh  size  to  be  L  *  0.578".  From 

Figure  9  it  is  found  that  the  total  length  of  one  thread  is  » 

1.06  L  -  0.612".  The  perimeter  of  the  wire  is  Pe  =  0.91".  The  total 

2 

cross-section  of  an  individual  wire  was  found  to  be  A  *  0.0507  in  . 


r* - b  -  0.37— 


s  T~  i 

O  3 _ \ 


R»  0.075“ 
NOMEX 


COATING 


FLOW  DIRECTION 


WIRE  CROSS-SECTION 
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Finally,  a  crude  correction  to  the  calculated  stagnation  temperature 
has  been  made  using  a  value  1000°R  smaller  than  the  one  plotted  in 
Figure  20.  This  corrected  agrees  more  closely  with  the  measured 
values  reported  in  Reference  22. 

The  energy  equation  is  applied  in  the  following  form: 

Vpc  ~  -  A  •  h(Tn  -  T)  +  A  .  a  •  I 
da  conv  0  rad 

x  eol(Fr  .  T4  F ,  .W4  -  T4)(F,  .  \(T4  -  T4)  ^ 

K  (p-spjj^  sp  (p-ea) J  ea  \(p-sp)2^  sp  J 

where 

Volume  V  2  •  L  *  A 

tot 

conv.  Area  A  *  2  •  L  ^  •  Pe 
conv  tot 

Rad.  Area  A  ^  "  2  *  L  •  h  (projected  area) 

The  shape  factors  are  taken  as  the  following  (without  further  justi¬ 
fication)  : 

outer  surface 


F(p-sp)1  ‘  2/3  J 


F  .  »  1/2  V  inner  surface 

IP-«P>1  J 


To  compare  bett«r  with  the  real  stagnation  temperature,  a  crude 
correction  has  been  supplied  by  using  a  value  10Q0°R  smaller  than  the 
one  shown  in  Figure  5. 

Figure  8  shows  the  plots  for  the  eight  cases 

^To  y0  -  1000  ^ 

lariu^  turb..  lam.  turb. 

/\.  /\.  /\.  /\. 

rad  rad  rad  rad  rad  rad  rad  rad 

Abbreviations:  T  •  turbulent 

L  =  laminar 

R+  »  radiation  increased 


lire  radiation  terms  seem  to  have  little  influence  for  t  <  10  sec 
and  for  the  decreased  stagnation  temperature,  where  the  change  in  the 
radiation  term  affects  T  only  by  less  than  10°R.  Comparison  with  the 
reported  measurements  indicates  that  only  a  more  exact  calculation,  i.e., 
including  conduction  effects  inside  the  wire,  will  give  results  in  the 
vicinity  of  the  available  date. 
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Figure  ltt. 


lo  ,i  i  :  Data  f or  Goodyear  .'rajertory 


T<0)  -  44.2 S  Me. 


3100 


T 


Figure  20. 


TABLE  VII 


Data  from  Flight  SP-5 


Launching  time  6:57:0.0 

Given  point:  at  121,000  ft  altitude  Mach  5.5 


Launching  Data 


time 

altitude 

V1 

al 

M1 

Pl^O 

P] 

[atm] 

[T  +  ... 
f  c  ] 

[1000  ft] 

[ft/sec] 

[ft/sec] 

44.28 

118 

5665 

5.458 

0.02576 

5.36  x  10"^ 

45.28 

123.1 

5585 

1063 

5.254 

0.02775 

4.05  x  10  ^ 

46.28 

128.2 

5516 

1088 

5.069 

0.02981 

3.30  x  10  ^ 

47.28 

133.3 

546w 

1112 

4 . 910 

0.03174 

2.77  x  10  * 

48.28 

138.4 

5405 

1136 

4.758 

0.03377 

2.32  x  10  * 

49.28 

143.5 

5360 

1161 

4.617 

0.03582 

1.95  x  10_^ 

50.28 

148.6 

5320 

1186 

4.486 

0.03791 

1.67  x  10 

51.28 

154 

5275 

1207 

4.370 

0.03991 

1.39  x  lO'^ 

52.28 

159 

5240 

1227 

4.271 

0.04175 

1 . 195x  10  * 

53.28 

164 

5205 

1238 

4.204 

0.04306 

1.03  x  10  J 

54.28 

169 

5170 

1244 

4.156 

0.04404 

0.89  x  10'.' 

55.28 

173 

5140 

1247 

4.122 

0.04476 

7.95  x  10 

56.28 

179 

5110 

1248 

4.095 

0.04534 

6.75  x  10  7 

57.28 

184 

5080 

1248 

4.07] 

0.04586 

5.9  x  10*? 

5.13  x  10'* 

58.28 

189 

5050 

1247 

4.050 

0.04633 

59.28 

193.9 

5019 

1243 

4.038 

0.04660 

4.45  x  10  . 

60.28 

198.8 

4989 

1238 

4.030 

0.04678 

3.85  x  10"? 

61.28 

203.5 

4958 

1232 

4.024 

0.04691 

3.38  x  10 

62.28 

208.4 

4928 

1222 

4.033 

0.04671 

2.91  x  10'* 

63.28 

213.4 

4897 

1209 

4.050 

0.04633 

2.50  x  10'* 

64.28 

218.1 

4867 

1196 

4. 0b9 

0.04590 

2.15  x  10 

65.28 

222.8 

4836 

1182 

4.091 

0.04542 

1.87  x  10~* 

66.28 

227.4 

4804 

1166 

4.120 

0.04480 

1.59  x  10  * 

67.28 

232.1 

4773 

1148 

4.158 

0.04400 

H4 

X 

O 

1 

"  £ 

68.28 

236.8 

4742 

1128 

4.204 

0.04306 

1.14  x  10 

69.28 

241.5 

4712 

1105 

4.264 

0.04188 

9.55  x  10  * 

8.00  x  10 

6.78  x  1<T* 

70.28 

246.0 

4681 

1077 

4.346 

0.04035 

71.28 

250.  3 

4650 

1047 

4.441 

0.03867 

72.28 

254.6 

4617 

4,540 

0.03703 

5.60  x  10  j? 

73.28 

2  56.9 

4584 

985 

4.660 

0.03518 

4.55  x  10~i 

74.28 

263.2 

4553 

965 

4 . 720 

0.03431 

4.70  x  10 

70 


TABLE  VI1.I 


Trajectory  Parameters 


time 

[sec  j 

M 

a2 

T„/T 

xJ  25 

44.28 

0.4094 

1.0335 

4r>.  28 

0.4119 

1.0339 

7l-.  28 

0.4142 

1.0343 

4  '  ,28 

0.4165 

1.0347 

48.28 

0.4189 

1.0351 

89.28 

0.4214 

1.0355 

50.28 

0.4238 

1.0359 

51.28 

0.4262 

1.0363 

52.28 

0.4283 

1.0367 

53.28 

0.4298 

1.0369 

54.28 

0.4310 

1.0372 

85.28 

0.4318 

1.0373 

;,6.23 

0.4325 

1.0374 

57.28 

0.4331 

1.0375 

58. 28 

0.4336 

1.0376 

59.28 

0.4340 

1.0377 

.0,28  | 

0.4342 

1.0377 

61.23 

0.4343 

1.0377 

62.28 

0.4341 

1.0377 

63.28 

0.4336 

1.0376 

64.28 

0.4331 

i  -375 

65.28 

0.4326 

^  0374 

66.28 

0.4319 

1.0373 

67.28 

0.4311 

1.0372 

68.28 

0.4298 

1.0370 

69.28 

0.4285 

1.0367 

70.28 

0.4267 

1.0364 

-1.28 

0.4247 

1.0361 

7  2.28 

0.4228 

1.0358 

73.28 

0.4206 

1.0354 

74.28  1 

0.4196 

1.0352 

VTi  to  po 

[°R!  [ lbf / ft.2  j 


6.734  3113.15  440.5 

6.307  3054.21  309.0 

5.935  3005.57  234.35 

5.626  2975.93  184.75 

5.340  2943.15  145.44 

5.083  2915.33  115.25 

4.851  2892.80  93.257 

4.651  2880.45  73.732 

4.484  2870.14  60.594 

4.373  2850.24  50.639 

4.295  2824.26  42.782 

4.240  2801.03  37.601 

4.197  2780.73  31.517 

4.159  2755.82  27.236 

4.125  2725.85  23.441 

4.107  2701.94  20.216 

4.094  2675.04  17.423 

4.068  2636.77  15.254 

4.099  2629.20  13.189 

4.125  2609.42  11.423 

4.155  2587.79  9  9162 

4.191  2560.75  8.7159 

4.237  2524.28  7.5134 

4.298  2486.00  6.4472 

4.373  2440.73  5.6047 

4.472  2378.43  4.8274 

4.610  2313.51  4.1973 

4.773  2225.48  3.7117 

4.946  2148.70  3.2015 

5.160  2062  89  2.7380 

5.270  1914.96  2.28 


7 1 


REFERENCES 


1.  Hritzay,  D.,  and  Wiant,  "Wire  Cloth  Structure  for  a  Radiating  Re- 
Entry  Vehicle,"  AVCO-Everett  Research  Report  123,  1962. 

2.  Kyser,  A.  C. ,  "The  Rotornet:  A  High-Performance  Hypersonic  Decelera-- 
tor  for  Planetary  Entry,"  NASA  CR-247,  June  1965. 

3.  Scott,  C.  J.,  and  Eckert,  E.  R.  G. ,  "Experiments  on  the  Thermal  Per¬ 
formance  of  Ribbon  Parachutes,"  Air  Force  Flight  Dynamics  Laboratory 
Report  AFFDL-TR-64-192,  May  1965. 

4.  Edwards  R.  C. ,  et.  aL,  "Study  and  Exploratory  Free-Flight  Investi¬ 
gation  of  Deployable  Aerodynamic  Decelerators  Operating  at  High  Alti¬ 
tudes  and  at  High  Mach  Numbers , "FDL-TDR- 04-35 ,  July  1964. 

5.  Nebiker,  F.  R. ,  "Aerodynamic  Deployable  Decelerator  Performance-Evalu¬ 
ation  Program,"  AFFDJ -TR-65-27,  August  1965. 

6.  Koh,  J.  C.  Y,,  and  Hartnett,  J.  P. ,  "Measured  Pressure  Distribution 

l| 

and  Local  Heat  Transfer  Rates  for  Flow  over  Concave  Hemispheres, 

ARS  Journal,  Vol.  31,  No.  1,  January  1961. 

7.  Korst,  K.  H. ,  Chow,  W.  L. ,  and  Zumwalt,  G.  W. ,  "Research  on  Transonic 
and  Supersonic  Flow  of  a  Real  Fluid  at  Abrupt  Increases  in  Cross 
Section  (Final  Report),"  University  of  Illinois  Eng.  Exp.  Station 

TR  392-5,  December  1959. 

8.  Cornish,  R.  H. ,  Ahimaz,  F.  J.,  Beadle,  C.  W. ,  and  Foster,  K. ,  "Mass 
Transfer  Cooling  of  Parachute  Materials,"  WADC-TR-58-684,  September  1959. 

9.  Hermann,  R. ,  "Problems  of  Hypersonic  Flight  at  t lie  Re-Entry  of  Satellite 
Vehicles,"  University  of  Minnesota  Rosemount  Aeronautical  Labora¬ 
tories  RAL-RR-153,  November  1958. 

10.  Eckert,  E,  R.  G. ,  Hartnett,  J.  F. ,  and  Tutuacliff ,  R.D.,  "Transient 
Temperatures  of  Parachutes  During  Descent,"  WADC--TN- 57-320 ,  August  1957. 

11,  Eckert,  E.  R.  G. ,  "Survey  on  Heat  Transfer  at  High  Speeds,"  WADC-TR-54- 
70,  April  1954. 

12.  Rose,  P. ,  Pro! stein,  K. ,  and  Adams,  M.  S. ,  "Turbulent  Heat  Transfer 
Through  a  Highly-Cooled,  Partially  Dissociated  Boundary  Layer,"  journ. 
of  Aero.  Sciences  Vol.  25,  No.  12,  December  1958. 

13,  Detru,  R.  W. ,  and  Hidalgo,  H.  ,  "Generalized  Heat  Transfer  Formulas 

a i-d  Graphs  for  Nuse  Cone  Reentry  into  Atmosphere,"  Journ.  American 
Rocket  Society,  Vol.  31,  1961. 


72 


14.  Eckert,  E.  R.  G. ,  et  al,  "An  Experimental  Study  of  the  Distribution 

of  Convective  Heat  Transfer  to  a  Large-Scale  Model  of  Parachute  Cloth, 
AFFDL-TR-66-13 ,  May  1966. 

15.  Goldstein,  S.,  Modern  Developments  in  Fluid  Dynamics,  Vol.  2, 

Oxford  University  Press,  N  w  York,  1938. 

16.  Sibulkin,  M.,  "Heat  Transfer  to  an  Incompressible  Turbulent  Boundary 
Layer  and  Estimation  of  Heat  Transfer  Coefficients  at  Supersonic 
Nozzle  Throats,"  J.  Ae.  Sciences,  Vol.  23,  No.  2,  February  1956. 

17.  Kudryavtsev,  Y.  V.,  et  al,  Unsteady  State  Heat  Transfer.  Iliffe 
Books  Ltd.,  London,  1966. 

18.  Schneider,  P.  J.,  Conduction  Heat  Transfer.  Addison-Wesley  Pub¬ 
lishing  Company,  Inc,,  Reading,  Mass.,  1957. 

19.  Carslaw,  H.  S.,  and  Jaeger,  J.  C.,  Conduction  of  Heat  in  Solids, 
Oxford  University  Press,  1950. 

20.  Schneider,  P.  J.,  Temperature  Response  Charts.  John  Wiley  and  Sons, 
Inc,,  New  York,  1963. 

21.  Heisler,  M.  P,,  "Temperature  Charts  for  Induction  and  Constant 
Temperature  Heating,"  Trans.  ASME,  Vol.  69,  April  1947. 

22.  Bloetscher,  F. ,  and  Arnold,  W.  V.,  "Aerodynamic  Deployable  Decelera- 
tor  Performance-Evaluation  Program,"  AFFDL-TR-67-60,  December  1967. 


23.  Noreen,  R.  A.,  and  Hosker,  R.  P.,  "Pressure  Distribution  on  Several 
Rigid  Hyperflo  Models  at  Mach  3.0  and  4.0,  AFFDL-TR-67-84,  June  1967. 


UNULAUo]  MED 


Sr i  unt v  CI.ISM fi<  i*t »<wi 


DOCUMENT  CONTROL  DATA  •  R  &  D 

{ Sec  nr  tty  <  la  v  s.1  fnati  tm  of  tit  ft-,  hmU  of  ahstra «  t  tit  nf  tt'shr  *  mg  ,»nn.>  httii'ti  mu*«?  hr  rnfern/  rt/irn  thn  y-r.itl  rpin-ff  J  ■»  i  *»///«♦  if} 


«  on  I  GlN  A  TING  AC  rrv'ITv  t' Torpor*  10  rti.l/ior)  MM'l'HT  ft  1  Ciiril  l»  (  l  H  K ft  I  »  I  I  A  III  IN 

University  of  Minnesota  j  Unclassified 

Institute  of  Technology  hr  lUTUTTT - " 

N/A 


Mechanical  Engineering  Dept,  „ 

jSh^s^s^Mhpsast^Sli^ _ _ _ - _ M _ 

3  REPORT  TITLE 

Tiie  Prediction  of  Material  Temperature  on  Woven  Retardation  Devices 


A  DESCRIPTIVE  NOTES  (Typo  o/  r?pof.'  oitu*  if*c/ii*ivp 

Final  Report  August  IS 66  -  September  196? 


5  Au  THORlS)  (Fir*  T  narntf,  middle  initial.  /*w.  name) 

Scott,  Charles  J 


...  — !'F35^r^c.102a 

b.  PROJEC  T  NO  6065 

c.  Task  No.  606505 


’0.  DtSTRI  OU  T  ION  STATEMENT 

This  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  the 
AF  Flight  Dynamics  Laboratory  (FDFR),  Wright-Patterson  AFB,  Ohio. 


12.  SPONSOR!:"".  MILITARY  ACTIVITY 


7a.  TOT  At.  ML.'.  OP'  P  A  G  IE  5 

73 

7b.  NO 

OF  M  r  F  * 

?3 

S«.  ORIGINATOR'S  REPORT  NUM0ERUI 

University  of  Minnesota 

HTL-TR-71 

nb.  OTHER  REPORT  NOIS)  (Any  other  numbers  that  may  be  an  signed 
this  report) 

AFFDL-TP.  -67-170 

AFFDL  (FDFR) 

Wright- Patterson  AFB,  Ohio 


13  ABSTRACT 


The  primary  purpose  of  this  report  is  to  present  analytical  techniques  for  the  pre¬ 
diction  of  instantaneous  parachute  temperatures  for  any  prescribed  flight  condition 
and  history  of  operation.  Sample  calculations  are  presented  which  illustrate  the 
technique  by  comparing  measured  parachute  temperature  histories  with  several  pre¬ 
dictive  techniques.  A  variety  of  exact,  transient  heat  conduction  solutions  are 
included  (for  a  cylindrical  geometry)  which  cover  a  range  of  trajectory  possibilities, 
This  was  done  because  the  transfer  of  heat  +~  -  ''cnr-rtional  fabric  or  coated-fabric 
parachute  is  governed  primarily  by  the  large  resistance  to  heat  conduction  into  the 
interior.  The  commonly  used  infinite  tnermal  conductivity  case  is  shown  '«.  >  prc..uce 
surface  temperatures  which  fall  considerably  below  the  surface  temperature,  predicted 
by  more  exact  (complete)  energy  balances.  For  the  extreme  applied  heat  flux  case 
(atmospheric  entry),  internal  temperatures  are  accurately  predicted  by  assuming  the 
wail  temperature  instantaneously  assumes  a  steady  state,  radiation  equilibrium  value 
(qw  *  0)  and  internal  conduction  governs  the  process.  An  approximate  method  is  pre¬ 
sented  for  the  calculation  of  parachute  fabric  temperatures  when  both  the  internal 
and  external  resistances  to  the  flow  of  heat  are  of  importance.  The  results  arc 
comfore  with  flight  tes  data  on  a  full-scale  parachute. 

The  distribution  of  this  Abstract  is  unlimited. 
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